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Abstract 

We give a determination of the equivalence group of Euler-Bernoulli equation 
and of one of its generalizations, and thus derive some symmetry properties of 
this equation. 
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1. Introduction 

Consider a collection T of differential equations of the form 

A(x,y (n) ;^)=0, (1) 

where x = (x , . . . , x p ) is the set of independent variables and ?/(„) denotes y and 
all its derivatives up to the order n, while the parameter A denotes collectively 
the set of all arbitrary functions specifying the family element in J 7 . On the 
other hand, let G be the Lie pseudo-group of point transformations of the form 

x = <p{z,w), y = if)(z,w), (2) 

where z = (z , . . . , z p ) is the new set of independent variables, and w = w(z) is 
the new dependent variable. As Tresse explained in [1, P. 11], the elements of 
G depend in general on arbitrary functions and not on arbitrary constants, and 

•Tel: +27 11 717 6237; Fax: +27 86 579 6248 
Email address: jean-claude .ndogmoOwits . ac . za (J.C. Ndogmo) 



Preprint submitted to Nonlinear Analysis: Real World Applications 



October 28, 2011 



G is therefore infinite-dimensional. We say that G is the group of equivalence 
transformations of (fT]) if it is the largest Lie pseudo-group of point transforma- 
tions that map ([T]) into an equation of the same form, that is, if in terms of the 
same function A appearing in (flj, the transformed equation has an expression 
of the form 

A(z,w (n) ;B) = 0, (3) 

where B denotes collectively the new set of arbitrary functions of the equation. 
Equivalence groups play an important role in, amongst others, the classification 
of differential equations by means of their invariant functions, and they are also 
a valuable tool for the identification of these invariant functions 

In this paper we obtain the equivalence transformations of Euler-Bcrnoulli 
equation (EB). The difficulty with this equation is due to the high order of 
derivatives that it involves, especially in the two-dimensional case, and this 
makes the standard method of finding the equivalence transformations by direct 
analysis intractable. We therefore recourse to the transformation of differential 
operators that we discuss in the next section. Next, we proceed to the determi- 
nation of the equivalence group of (EB) and of one of its generalizations, and 
derive some symmetry properties of this equation. This includes the full sym- 
metry group of the equation, which was obtained only partially in [6| in terms 
of some undetermined set of functions. 



2. Transformation of differential operators 

An expression of differential operators in terms of new variables yields an 
effective algorithm for implementing the transformation of a differential equation 
under a given change of variables, and this can also greatly simplify calculations, 
based on the properties of these differential operators. 

Suppose for instance that the general change of variables ([2]) that maps |lj 



2 



into (J3|) is given more explicitly in the form 



x l = tpi(z,w), z — (z , . . . , z p ), w = w{z) (4a) 
u = tp{z,w) = T(z). (4b) 



Then the last equality gives 

a f = X^"1TT' fori = l,...,p, (5) 
az 1 ^— ' 02 l 

and for invertible transformations the linear system ([5]) yields solutions of 
the form 

^ T ,dT(z) , . 1 
"< 2J V ":— ]77-- forj = l,...,p, 

i=i 

for some functions ^ — ^j(z). In other words, an expression of differential 
operators in terms of the new independent variables is given by 

'), >I' ; . where ^ % il ■ ( 6 ) 

t 

We shall make use of some properties of the linear differential operators ^ to 
simplify calculations, especially when dealing with high-order derivatives. 

3. Determination of the equivalence group 

We now move on to consider in this section the problem of finding the 
equivalence group for the one-dimensional Euler-Bernoulli equation, which is a 
model for calculating the load-carrying and deflection characteristics of beams. 
Since the late 19th century, following the successful demonstration of this theory 
in the constructions of the Eiffel Tower and the Ferris wheel, Euler-Bernoulli 
equation (also known as Engineer's beam theory) became the cornerstone of 
engineering. 

This equation can be put in the form 

u 

/to-H) + m (z) J k2=°> <>0 > 0<x<V, (7) 



dx 2 V dx 2 I dt 2 
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where f(x) > is the flexural rigidity, m(x) > is the lineal mass density and 
u(t, x) is the transversal displacement at time t and position x from one end of 
the beam taken as origin. 

In the relatively recent research literature, Euler-Bernoulli equation, which 
is often simply referred to as Euler-Bernoulli beam or beam equation, has been 
frequently discussed. In particular, Gotllieb 7] has investigated the iso-spectral 
properties of this equation and its non-homogeneous variants in connection with 
the unit beam (i.e. / = 1 and m — 1). More recently, Wafo Soh 6] considered 
the equivalence problem for Euler-Bernoulli beam from the Lie symmetry ap- 
proach, while Morozov and Wafo Soh [8| investigated the same problem using 
Cartan's equivalence method. 

It is clear that in any change of variables that preserves the form of (J7J), the 
transformation of x should depend on a single variable (due to the form of / 
and m), while the linearity of the equation forces the transformation of u to 
be linear, and that for the independent variables not to involve the dependent 
variable. We are therefore led to look for equivalence transformations of (0 in 
the form 

t = R(y,z), x = S(z), u = L(y, z)w + J(y, z) (8) 

Under (|SJ), the coefficient 7 i of w yyyy in the transformed equation takes the 
form 

7l = (fB*Ly(B*Si) 

and its vanishing implies that R = R(y). Inserting this new expression for R 
into ([8} , the vanishing of the coefficient of w y in the new transformed equation 
yields the condition 

"2RyLy LRyy = 0, 

and thus L = h(z)-JR^, for some function h. If we now transform ([7|) using this 
additional information and write down the coefficient 72 of w z as a polynomial 
in f(S) and its derivatives, then the coefficient <5i of fs,s in this expansion is 

5 l = S 4 z (2h z S z - hS zz ), 
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and its vanishing gives h = fci%/5z", where k\ is a constant, so that u = 
h^RySzf^w + J. With this expression for h, the vanishing of the coefficient 
of fs in 72 leads to the condition 



^ S zz ~l~ 2 S z S Z zz — 0, 

and this in turn shows that 

x = S = (k 2 z + k- 3 )/(k 4 z + 1), k 2 — k 3 k 4 ^ 0, 

for some constants k 2 , . . . ,k 4 . With all the information obtained up to this point 
on R, S, and L, the coefficient 73 of w now takes the form 

4(fc 2 - fc 3 fc 4 ) 8 m (3flg„ - 2R y R yyy ) 
73 ~ (1 + k 4 z) ie 

showing that R must satisfy an equation similar to that for S 1 , namely 

Consequently, we must have 

t = R=(k 5 y + k 6 )/(k 7 y+l), k 5 ~k 6 k 7 ^0, (9) 

for some arbitrary constants A; 5, . . . , k 7 . With this new information and all others 

obtained thus far, ([8]) now takes the form 

k 5 y + k 6 k 2 z + k 3 

t=— — , x = — (10a) 

k 7 z + 1 k 4 z + 1 v ' 

, (fc 2 -fc 3 fc 4 ) 1 / 2 (fc 5 -fc 6 fc 7 ) 1 / 2 

u = h — (Ttwtm — w + J(2/ ' z) - (10b) 

Under (fT0|) . the transformation of Euler-Bernoulli beam retains a constant 
component, that is, a component (or sum of terms) free from w and its deriva- 
tives, and which appears to contain the only terms involving J and its deriva- 
tives. In fact this component is simply the transformation of u — J under the 
change of independent variables (|10a[) . and it is readily found that the vanishing 
of this constant component amounts to setting 

fc 8 - k 10 k^y + k 4 (-kg + kgkfy) + k\ (-k 9 z + k 14 k^yz) 
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where fcg, . . . fcn are some additional constants of integration. 

If we now perform a transformation of ([7J under the change of variables (|10[) 
with J as in (|11[) . the transformed equation takes on the form 



<9 2 / d 2 w 



where 



F(y,z) 
M(y,z) 



dz 2 



(1 + k 4 zf 



M(y,z\ 



d 2 w 
dy 2 



(12a) 



{l + k w ){l + k iZ ) S{S)l 
{l + k 7 y) 3 (l + k i z) 3 (k 2 - 



■ m(S). 



(12b) 
(12c) 



(fc s -fe 6 fe 7 ) 2 (l + fc^) 7 

Therefore, the change of variables given by (ITU1) and (|TT|) will represent an 
equivalence transformation of ([7]) if and only if the functions F and M in (|12[) 
do not depend on y, and this is possible if and only if k-j — 0. However, as in (|11[) 
the value kj — is not allowed, we must recompute the appropriate expression 
for J when k-j = 0. So, letting kj = in ^ and finding the corresponding value 
for J gives 

-fc - k 9 y + fc 4 (fc 8 + fcioy) + kjih + k w y)z 



J 



(13) 



fc4(l + /S4Z) 

for some new constants ko , and k% , . . . , fcio , and the transformed equation now 
takes on the required form 



d 2 



d 2 w 
dz 2 



M{z) 



d 2 w 
dy 2 



0, 



(14a) 



where 



F(z) = {l + k iZ ) 1 



fc|(l +/c 4 z) 3 



(fcg - hk 4 )k 5 
(1 + A: 4 z) 2 

(fc 2 - k 3 k A )k 5 
(1 + fc 4 z) 2 



1/2 



/(£), 



1/2 



(14b) 
(14c) 



We have thus obtained the following result. 
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Theorem 1. The group G of equivalence transformations of Euler- Bernoulli 
equation is given by the equivalence transformations 

t = k 5 y + k 6 , x= k f Z + ks (15a) 
kiz + 1 

f(k 2 - k 3 k±)k 5 \ 1/2 
u = ki [ n | - ] w + J, (15b) 

where 



J = 



, (l + k 4 z) 2 , 
-ko - k g y + k 4 (k s + k w y) + fcf (fc 8 + k w y)z 



ki(l + k^z) 

is the function obtained in (113|) . and the transformed equation is given by (|14[) . 



We also note that in the successive transformations of (|7|) which led to (|T2|) , 
the arguments of the arbitrary function / and m were treated as dummy, by 
being ignored. Thus, in view of (|12[) . we readily obtain the following result for 
an extended version of Euler-Bernoulli Equation. 

Lemma 1. For the generalized version of Euler-Bernoulli equation of the form 
(>1 f f(t, x )^)+ m (t, x )^=0, (16) 



dx 2 \ ' dx 2 J ' dt 2 

in which f and m are assumed to be functions of both t and x, an equivalence 
subgroup is given by the equivalence transformations (|10j) and (|llj) . and the 
transformed equation takes on the form (|12[) . 

We may not however assert as yet that (fT0| and (fTTj) define the equivalence 
group of (|16j) . because due to the form of / and m in JJJ, these transformations 
were obtained under the assumption that x = S(z), depends only on a single 
variable. But this restriction in no longer a priori permitted in the case of (|16p 
where / and m are functions of both t and x. 

The equivalence group of (fTl))) should thus be looked for in the form 

t = R(y,z), x — S(y,z), and u = L(y, z)w + J(x, y), (17) 

where as usual, the functions R, S, L and J are to be found. It first follows from 
the lemma and the invertibility of p7|) . that 

R y S z L ^ 0, wL ^ 0, (18a) 
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where 



W — RyS Z — R Z Sy. 



(18b) 



It is also clear that in the transformed equation, terms of the highest order 
may only come from the transformation of u xxxx . Moreover, to find such terms, 
we may assume that u — Lw in (|17[) , and since L ^ 0, by Leibnitz rule we may 
assume for such purpose that u = w(y, z). Since the term (R z l w ^) w yyyy appears 
in the polynomial expansion of the transformation of u XXXXl the vanishing of 
this term yields R = R(y). 

On the other hand, it follows from the formulas (j6|) that 

d t = -(S Z dy-Syd Z ), d X = -(-R Z d V +Ryd Z ), (l9) 

and for R = R(y), this reduces to 



1 „ S„ „ 1 



d t = — d y --^d z , d x = —5 z . (20) 



IXy IXyOz Oz 



Consequently, terms involving derivatives of w w.r.t. y in the transformed equa- 
tion can only come from the transformations of derivatives of u w.r.t. t, i.e. from 
the transformation of uu- So, transforming utt by setting u — L(y,z)w in (117|) 



shows that the term in w yz has coefficient 



-(2SyL)/(R 2 y S z ). 

The vanishing of this term gives S — S(z), and hence the transformations of 
(TTO|) must also be sought in the form 

t = R(y), x=S(z), u = L(y,z)u + J(y,z), (21) 

just as in the case of ([7]). Given that in the transformations of ([7} starting 
with a change of variables of the form (|2"Tj) . and leading to (|10l) . (jTTJ) , and (jl^j) . 
the arguments of the arbitrary functions / and m are ignored, and because ([7J 
and (1161) differ only by these arguments, we may conclude that the equivalence 
transformations of (fTTpj) must also be of the form dTU|) - (TTTj) . Since the resulting 
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transformations of (|16p which are also given by (fT2")l are of the required form, 
we have thus obtained the following result. 

Theorem 2. The equivalence transformations 

k 5 y + k 6 k 2 z + k 3 

t = — — , x — — — 

k-jz + 1 K4Z + 1 

, (fc 2 -fc 3 fc4) 1/2 (fc 5 -fc 6 fc 7 ) 1 / 2 

u=k i T\ 1 u \7T~TT. \ w + J{y,z) 

(1 + k 7 y)(l + k 4 z) 



where 



^ _ k s - k w k%y + fc 4 (-fc 9 + kntfy) + k\ (-k 9 z + kgkfyz) 
kikj(l + kjy)(l + k^z) 

and which are given by (|TU)) and pT|) . define the equivalence group of the 
extended form (|16|) 0/ ([7]) . 

4. Symmetry properties 

It is clear that the equivalence group of a differential equation contains ac- 
cording to our definition the largest symmetry group G e of the equation that 
is free from the arbitrary functions of the equation. However, such a symmetry 
group is only a subgroup of the so-called principal Lie group, which is the largest 
symmetry group that holds for all arbitrary functions, but which may depend 
on these arbitrary functions. 

Symmetries of Eulcr-Bcrnoulli equation were calculated in [6| , but were ob- 
tained only in a conditional form, and under the assumption that the functions 
/ and m satisfy certain unsolved complicated equations. 

It appears from the form of (TI4]) that (TT5)) will represent a symmetry trans- 
formation of (0 only if the transformation of x is the identity transformation, 
x — z. This means that we have to assume k± = in (|15p. and because the 
corresponding expression for the denominator of J vanishes for ki = 0, we have 
to recompute the appropriate value of J. To do so, we transform ([7J using (|10[) . 
with kj = &4 = &3 = 0, and ki — 1. It then follows that ([7]) will be invariant 
under these transformations if and only if fc 5 = 1 , and 

J = VP2 + Pa + z (ypi + p 3 ) . (22) 
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We thus have the following result. 

Theorem 3. The largest symmetry group of Euler- Bernoulli equation that is 
independent of its arbitrary functions is a six-parameter group G e given by the 
transformations 

t^y+pQ, x = z, and u = p$w + J, (23) 
where J is given by (|22[) , and with corresponding symmetry generator 

v =Ped t + (P4+P2t + +P3X + pxtx +p 5 u) d u , (24) 
where pi, . . . ,pq are the group parameters. 

We note that the function J in (|2"2"]) is the fundamental solution of Euler- 
Bernoulli equation that does not depend on its arbitrary functions. Owing to 
the linearity of (O, if we replace the function J in (|2"3")) by any other solution 
of ((7|), the resulting transformation remains a symmetry Q, [hJ], but which 
may however depend on the arbitrary functions of the equation, and this is 
contrary to our assumptions, and for practical considerations such solutions are 
not always available. 

On the other hand, due to the superposition principle of linear equations, 
(|23p with J replaced by an arbitrary solution S of ([7|l. remains a symmetry of 
Euler-Bernoulli equation, and one readily verifies that this is the most general 
symmetry group of this equation. 
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